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The geometric, kinetic, and total inductances and the attenuation constant are theoretically ana- 
lyzed for a thin-film superconducting coplanar waveguide (CPW) resonator consisting of a current- 
carrying central conductor, adjacent slots, and ground planes that return the current. The analysis 
focuses on films of thickness d obeying d < 2 A (A is the London penetration depth), for which the 
material properties are characterized by the two-dimensional screening length A = 2A 2 /d. Introduc- 
ing a cut-off procedure that guarantees that the magnitudes of the currents in the central conductor 
and the ground planes are equal, new and simpler results are obtained for the kinetic inductance and 
the attenuation constant for small A. Exact results for arbitrary A are presented for the geometric, 
kinetic, and total inductances in the limit of tiny slot widths, and approximate results are presented 
for arbitrary slot widths. 
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I. INTRODUCTION 

Superconducting coplanar waveguide (CPW) res- 
onators have found increasing use in microwave inte- 
grated circuits. 1 Their applications include microwave 
filters in mobile and satellite communication, 5-9 
superconducting qubits, 10-14 circuit quantum 
electrodynamics, 15,16 radiation detectors, 17-22 and 
parametric amplifiers. 24 Important characteristics of 
a CPW resonators are its inductances (geometric, 
kinetic, and total) and attenuation constant. Although 
well-known expressions for these characteristics are in 
common use, 1-3 it is the purpose of this paper to review 
the assumptions underlying their derivations and to 
derive new expressions when these assumptions cannot 
be justified. 

Considered here is a CPW resonator fabricated from a 
thin superconducting film of thickness d centered on the 
xz plane. The central conductor (\x\ < a, \y\ < d/2) is of 
width 2a centered on the z axis, and the superconducting 
ground plane occupies the regions \x\ > b, \y\ < d/2. 
Slots lie on either side of the central conductor in the 
gaps at a < \x\ < b. 

Section II reviews the derivation of well-known results 
for the geometric, kinetic, and total inductances at low 
frequencies when the current is essentially all supercur- 
rent. The current density in the central conductor and 
ground planes is assumed to be far below the depairing 
current density, so that the current dependence 25 of A is 
negligible. When all dimensions, including the thickness 
d, are much larger than the London penetration depth A, 
there is no penetration of the self-field into the conductors 
except within a distance A from the surface. However, 
when d < 2 A, the screening of the self- field is no longer 
governed by A but by the Pearl length 26 A = 2X 2 /d. 
In this case, perpendicular magnetic fields can penetrate 
into the film to a distance of the order of A from the film 
edges. When these depths of penetration (A or A) are 



much smaller than a, the method of complex magnetic 
fields 27 can be used to approximate the magnetic-field 
and supercurrent distributions. I use this approach to re- 
cover the well-known result for the geometric inductance 
and to derive a new, simpler expression for the kinetic 
inductance. Needlessly complicated expressions for the 
kinetic inductance were obtained in earlier derivations 
assuming cutoff procedures that implicitly violated the 
requirement that the magnitudes of the current in the 
central conductor and the return current in the ground 
planes are the same. 

Section III contains exact solutions for the magnetic- 
field and supercurrent distributions and the geometric, 
kinetic, and total inductances in an idealized CPW for 
which b — ¥ a and the screening is characterized by the 
2D screening length A. In Sec. IV, I consider the re- 
alistic case of b > a, introduce a form for the current 
distribution in the central conductor that approximates 
its behavior as A/a varies from zero to oo, and obtain 
approximate expressions for the geometric, kinetic, and 
total inductances. 

In Sec. V, I account for normal-fluid losses to derive a 
new, simpler expression for the CPW's attenuation con- 
stant. Section VI contains a brief summary and discus- 
sion of the results. 



II. INDUCTANCES FOR A THIN-FILM CPW 
ASSUMING NO FIELD PENETRATION 

In all applications of CPWs the central conductor car- 
ries an alternating current I z e lu}t . Except at very high 
frequencies or at temperatures close to the transition 
temperature T c , the current is essentially all supercur- 
rent, and the normal- fluid contribution due to the flow 
of thermally excited quasiparticles is negligibly small (see 
Sec. V). To avoid including the factor e lult for all quanti- 
ties proportional to the dynamic current, let us consider 
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the behavior at a time t for which e lut = 1. The central 
conductor (|x| < a, \y\ < d/2) carries a current I z in the 
z direction, and the semi- infinite ground planes (|x| > b) 
carry the return current in equal amounts. When the 
film thickness obeys d -C a and X 2 /d -C a, magnetic 
flux per unit length $' passes up through the right-hand 
slot (a < x < b) and down through the left-hand slot 
(—6 < x < —a). We assume that b — a >• A 2 /d, so that 
we can neglect the magnetic flux that penetrates into the 
film edges adjacent to the slots. See Figs. 1 and 2. 

The complex magnetic field 27 W{C) = H y (x,y) + 
iH x (x, y) describing this situation is, with ( = x + iy, 



no 



%A 



(C 2 -a 2 )V2((2 -b 2 ) 1 ^' 



(1) 



where A is a constant determined below. The magnetic 
flux density in the plane y = in the slots (a < \x\ < b) 
is B y (x,0) = fioH y (x,0), where 



Hy(X,0) = 



Ax 



x\y/ \x 2 - a 2 ){b 2 - x 2 )' 



(2) 



and the average sheet-current density in the central con- 
ductor and the ground plane, averaged over their thick- 
ness, is given by K z (x) = —2H x (x, e), where 



K z (x) 



2A 



^(a 2 -x 2 )(b 2 -x 2 ] 
2A 



^/(x 2 -a 2 )(x 2 -b 2 Y 



, \x\ < a, (3) 
'x\ > b. (4) 



[Figure 1 shows plots of H y (x,0) and K z {x) calculated 
from Eqs. (2)-(4) for a = 1, b = 2 and A = 1.] For 
the general case, the constant A can be determined from 
the requirement that the current carried by the central 
conductor is I z , which yields the relation 



(4A/b)K(k), 



(5) 



where K(k) is the complete elliptic integral of the first 
kind of modulus k — a/b. The magnetic flux per unit 
length <£>' can be obtained by integrating Eq. (2) over 
the slot, and the result is 



Mo A 



K(k') 



^oK(k') 
AK(k) ' 



(6) 



where k' = vl — k 2 is the complementary modulus. 

The complex potential from which %{C) can be derived 
via H(() = dg(()/d( is 



G(C) = (iA/&)F(arcsin(C/a),a/6), 



(7) 



where F(<f>, k) is the elliptic integral of the first kind of 
amplitude <j) and modulus k. The contours in a contour 
plot of the real part of Q correspond to magnetic field 
lines, as shown in Fig. 2. These contours also correspond 
to contours of the vector potential A = zA z (x,y) (B = 
V x A), because A z (x,y) = —[ioRcQ(x + iy). 
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FIG. 1: H y (x,0) (solid) and K z (x) (dashed) vs x for a = 1, 
b = 2, and A = 1. 




FIG. 2: Contour plot of the real part of Q(x + iy) [Eq. (7)] 
as a function of x and y, where b — 2a when \ 2 /d <C a. The 
central conductor (|a;| < a) carries current I z and the ground 
planes (|x| > b) carry the return current equally. Contours 
correspond to magnetic field lines, which flow in a counter- 
clockwise sense around the central conductor. 



A. Geometric inductance for A 2 /d -C a 

From Eqs. (5) and (6) we recover the well-known result 
for the geometric inductance per unit length 28,29 



_ *' _ ^oK(k') 



h 4K(fc) ' 
which is well approximated by 

L' m ^^-[\n(4b/a)-(a/2b) 2 } 
for large b/a and by 



(8) 



(9) 



(10) 



3 




FIG. 3: L' m /fi , the geometric inductance per unit length in 
units of no, vs b/a, calculated for the case that A 2 /d <§C a from 
the exact expression Eq. (8) (solid) and the approximations 
Eq. (9) (dashed) and Eq. (10) (dotted). 



for (b — a) <C a. Shown in Fig. 3 are plots of L' m / [i 
calculated from Eqs. (8) (solid), (9) (dashed), and (10) 
(dotted). 



B. Kinetic inductance for \ 2 /d <C a 

For the calculations of this paper the current density 
is assumed to be everywhere far below the depairing crit- 
ical current density, such that only the weak-field (low- 
current) penetration depth A is involved. For this case, 
since the kinetic energy per unit volume of the supercur- 
rent is 30 no\ 2 j 2 /2, where j is the supercurrent density, 
the kinetic inductance per unit length L' k can be calcu- 
lated from 



1 x2 



dS, 



(11) 



where the integral is over the cross section of the CPW. 

When \ 2 /d <C a, the supercurrent density is given to 
good approximation by 

j z (x, y) = [K z (x)/2\] coshfo/A)/ sinh(d/2A) (12) 

except very close to the slots. Here K z (x) is given by 
Eqs. (3) and (4), and the integral of j 2 (x,y) over the 
film thickness is 



r-d/2 

J-d/2 



j 2 z (x,y)dy = K 2 z (x)q(d/X)/X, 



where 



q(u) — (sinhw + u)/8sinh 2 (u/2). 



(13) 



(14) 



Note that q(u) « 1/u when m«1, and q(u) s=y 1/4 when 
u > 1. 

Because K 2 (x) [Eqs. (3) and (4)] diverges at the edges 
of the slots, the integral of this quantity over x can be 



evaluated only approximately by cutting off the integrals 
at x = ±(a — 5 a ) and x = ±(b + 5b), where 5 a <C a 
and 5b <C b. The correct cutoff procedure is to choose 5 a 
and 5b such that the magnitudes of the current in the cen- 
tral conductor and the total return current in the ground 
planes are the same. Using this criterion to evaluate the 
currents for \x\ < a — 5 a and |a;| > b + 5b using Eqs. (3) 
and (4) leads to the condition that 

K-K^-M-^K) < 15 » 

or 5 a /a = 5b /b = e <C 1. Applying this cutoff procedure 
to evaluate the integral of K 2 {x) over x yields, to lowest 
order in e, 



L' u = 



li Xq(d/X) 



In 



4a(l- k)K 2 {k) Le(l + fc) 



2(1 -k) 



(16) 



where k = a/b. Equation (16) also can be written as 



L' k = ^q(d/X)g e (k,e), 



where 



9e(k,e) = 



In 



2(1 - k)K 2 (k) Le(l + fc) 



2(1 -k) 



(17) 



(18) 



is dimcnsionless and for typical values of a and b is of 
order unity (to logarithmic accuracy). 

In Rcfs. 29 and 31 the two cutoff lengths <5 a and 5b were 
assumed to be the same. This assumption violates the re- 
quirement that the calculated central-conductor current 
be the same as the magnitude of the return current in the 
ground plane. Both these papers 29 ' 31 derived expressions 
for L' k that contain two logarithmic terms arising from 
the cutoffs near a and b. Each of the two terms resembles 
the logarithmic term on the right-hand side of Eq. (16). 
Despite differences when b > a, the dependence of Eq. 
(16) upon a, b, k = a/b, and the cutoff parameter when 
(b — a) <C a is in approximate agreement with the more 
complicated Eq. (8) in Ref. 31 for the internal inductance 
of a superconducting coplanar waveguide. However, what 
we here call the cutoff lengths 5 a and 5b was called the 
stopping distance in Ref. 31. 

Watanabe et al. 29 did a similar calculation using w as 
the width of the center conductor, s as the slot width, 
and 5 = d/A as the single cutoff length. With the re- 
placements a — w/2, b — w/2+s, and 5 — d/A, the above 
expression for g f [Eq. (18)] agrees with the more compli- 
cated expression for g(s, w, d) given in Refs. 29 and 32 
when s <C w, except for a typographical error there: the 
k 2 in the denominator of the prefactor in their g(s,w,d) 
should have been k' 2 = 1 — k 2 . 

Mcservey and Tedrow 33 carried out a calculation of the 
kinetic inductance of a long, thin superconducting strip, 
and to handle a similar divergence they introduced a sin- 
gle cutoff length 5 (similar to 5 a or 5b) given by 5 w A/ 2 
when d > 2A and 5 ~ d/4 when d < 2 A. However, it is 
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now well known that for d < 2X, the equations governing 
the current density and self-field no longer contain A but 
the two-dimensional screening length (or Pearl length 26 ) 
A = 2X 2 /d. (See also Sec. V.) As a consequence, when 
d < 2 A, it is incorrect to assume 6 ~ d/4, and instead 
one should use cutoff lengths of order X 2 /d. 

As seen from Eqs. (8)-(10) and Fig. 3, for typical values 
of a and b, the geometric inductance per unit length L' m is 
of the order of hq. By comparison, the kinetic inductance 
per unit length L' k given in Eq. (17) depends strongly 
upon the ratio of d/X. When A <C d <C a, which satis- 
fies the assumptions leading to the results of Sec. II A, 
q(d/X) « 1/4 [Eq. (14)], and L' k « (p X/8a)g £ (k, e) « 
L' m . On the other hand, when A » d, q(d/X) w X/d, so 
that Eq. (17) becomes L' k w (p^A/ Aa)g c (k, e), in basic 
agreement with Eq. (3) of Ref. 29 when s -C w, except 
that the cutoff length should be S w A 2 /d instead of d/4, 
as assumed in Ref. 29. 



C. Inductances for arbitrary A 2 /d relative to a 

In this paper we explore how to determine both the 
geometric inductance and the kinetic inductance in thin 
films (d < 2A) for arbitrary values of A = 2X 2 /d rela- 
tive to the width 2a of the center conductor. In partic- 
ular, we are also interested in examining the condition 
A 3> 2a, when the kinetic inductance is expected to be 
much larger than the geometric inductance. Note that 
the equation L' k w (paA/4a)g e (k, e) cannot be used for 
this case, because it is a good approximation only when 
both d < A and e w A 2 /da < 1. Moreover, Eq. (8) for 
L' m is no longer valid, because its derivation required the 
assumption that A 2 /c? <C (6 — a), so that in the plane 
of the film the current-generated magnetic field was con- 
fined entirely to the slots. Therefore, to deal with thin 
films for arbitrary values of A = 2X 2 /d, including the 
case for which A > 2a, a new theoretical approach must 
be employed to obtain good estimates of L' m and L' k . 



III. THIN-FILM CPW WITH NARROW SLITS 

When the thickness d of a superconducting film cen- 
tered on the xz plane is somewhat less than the Lon- 
don penetration depth A and the film locally carries 
a supercurrent density j(x,y, z)e lut , solutions of the 
London equation 34 and Maxwell's equations reveal that 
j(x, y, z)e %ljjt is practically uniform across the thickness 
(i.e., independent of y). One then may deal only with 
the sheet-current density K = jd, so that the Lon- 
don equation becomes K = — (2/^ A)[A + (4>o/2tt)\7^], 
where A = 2X 2 /d, A is the vector potential in the 
superconducting film, (f>o = h/2e is the superconduct- 
ing flux quantum, and 7 is the phase of the supercon- 
ducting order parameter. 30 The gauge-invariant quantity 
A s = A + (0o/27r)V7, which is proportional to the su- 
perfluid velocity, is sometimes called the gauge-invariant 



vector potential. In Sees. Ill A and IIIB, I use this thin- 
film London equation and Maxwell's equations to obtain 
exact analytic solutions for the geometric and kinetic in- 
ductances per unit length of a CPW resonator with very 
small slit widths. As in Sec. II, to avoid including the 
factor e %u)t for all quantities proportional to the dynamic 
current, it is convenient to consider the behavior at a 
time t for which e tut — 1. 



A. Exact solution for a film with one long slit 

Consider a thin superconducting film in the xz plane 
with a long and narrow slit along the z axis containing N' 
flux quanta per unit length, such that the magnetic flux 
per unit length up through the film is $' = N'(f> n . We 
choose a gauge such that the vector potential has only a 
component along the z direction given by 



M{x,y) 



a(k) sin kx e~ k ^dk. 



(19) 



Similar expressions can be obtained for the correspond- 
ing magnetic flux density components Bo x (x,y) — 
dA (x,y)/dy and B 0y (x,y) — —dAo(x,y)/dx. The su- 
percurrent sheet-current density has only a component 
along the z direction given by the thin-film London equa- 
tion, 



K (x) = -(2/ Mo A)L4 (x,0) + (<W2tt) 7z 



(20) 



where J z (x) = dj(x,z)/dz — (7r$'/0o)sgn(x) with 
sgn(x) = 1 for x > and -1 for x < 0. According 
to Ampere's law, the sheet-current density is given by 
the discontinuity in B x across the film thickness. Using 
K (x) = [B x (x,0~) - B x (x,0+)]/no and Eq. (19), sub- 
stituting into Eq. (20), and taking the derivative with 
respect to x yields the equation 



f 

Jo 



a{k)k(l + k A) cos kxdk = -$'5(x), (21) 



where 5(x) is the Dirac delta function. Multiplying by 
cos k'x and integrating over all x yields 



a(fe) = - 



7rfe(l + fcA)' 



(22) 



The magnetic flux per unit length $' carried up 
through the film can be described briefly as follows: 
Bo y (x,0) is an even function of x/A that diverges at 
x = but decreases to zero as \x\ — >• 00. For y > (y < 0) 
and p = yjx 2 + y 2 > A, B a w p^'/np (B Q w -pQ'/np), 
where p = (xx + yy) / p. In other words, for y > 0, screen- 
ing by the superconducting film causes the magnetic flux 
density at large distances to look like that produced by 
a magnetic line charge. 

The vector potential has the following properties: 
Along the x axis, Aq(x,0) is an odd function of x/A, 
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FIG. 5: Sheet-current density K z (x) from Eqs. (27)-(29)], 
normalized to I z /a [Eq. (32)], vs x/a for A/a — 0.01 (solid), 
0.1 (dot-dashed), 1 (dashed), and 10 (dotted). 



FIG. 4: Contour plot of the vector potential A z (x,y) as a 
function of x and y, calculated from Eqs. (19) and (25) for 
the case of two narrow slits at x — ±a for A = a. The central 
conductor (\x\ < a) carries current I z and the ground plane 
(|x| > a) carries the return current. Contours correspond to 
magnetic field lines, which flow in a counterclockwise sense 
around the central conductor. 



and 4)(0,0) = 0. For p 



V^ 2 



y 2 > A, A a 



-(^'/^tan-^x/lyl), so that A (±oo,0) = =f$'/2. 
The corresponding sheet-current density is 



K (x) = - 



2$'sgn(x) 
npoA 



'(£)■ 



where 



/(«) - 



sin t dt 



t 



/(0) = tt/2, 4f(u)/d« 
for !t>l. 



-co at u 



t 2 + l ' 
0, and /(«) 



(23) 

(24) 
: 1/u 



B. Exact solution for a CPW with two narrow slits 

Consider now a superconducting film in the xz plane 
with two long and narrow slits parallel to the z axis, 
one at x = a containing N' positive flux quanta per unit 
length, and one at x = — a containing an equal number of 
negative flux quanta per unit length. The resulting vector 
potential, magnetic flux density, and sheet-current den- 
sity are obtained by a linear superposition of the reults 
obtained in Sec. Ill A: 



A z {x,y) 


= A (x - a,y) - A (x + a,y), 


(25) 


B{x,y) 


= B (x- a,y)- B (x + a,y), 


(26) 


K z (x) 


= K (x — a) — K (x + a) , 


(27) 




= (2/p A)[& -A z (x,0)\, \x\<a, 


(28) 




= -(2/p Q A)A z (x,0), \x\ > a. 


(29) 



These fields can be described as follows: B y (x, 0) is 
an odd function of x that diverges at x = ±a but de- 
creases to zero as \x\ —> oo. For p = \J x 2 + y 2 3> A 
and p » a, B y (x,y) w 4$'ax|y|/7r(a; 2 + y 2 ) 2 , and 
B x (x, y) w ±2$'a(x 2 - y 2 )/Tr(x 2 + y 2 ) 2 , where the upper 
(lower) sign holds for y > (y < 0). In other words, for 
y > 0, screening by the superconducting film causes the 
magnetic flux density at large distances to look like that 
produced by a dipole of magnetic line charges. 

The vector potential (see Fig. 4) has the following 
properties: Along the x axis, A z (x, 0) is an even function 

-(^'/^{tan- 1 ^- 
sjx 2 + y 2 > A 



of x. For p = yjx 1 + y 2 > A, A z 
a)/\y\]— tan~ 1 [(.T+a)/|?/|}, so that for p ■■ 
and p > a, A z (x, y) w 2^'a\y\/ir{x 2 + y 2 ). 

The corresponding sheet-current density K z (x) is an 
even function of x with discontinuities at x = ±o. For 

-4^'a/npox 2 ■ In the limit 



y 2 > A, K z (x) 



as A — > 0, K z {x) reduces to the form given by Eqs. (3)- (4) 
with b — a. See Fig. 5. 

i. Geometric inductance 

The geometric inductance per unit length L' m is defined 
via 



^ /■oo ^oo ^2 

^^m^ 2 = / / TT-dxdy. 

1 J-OO J-OG 



(30) 



-oo J — oo 



The two-dimensional integral over all space can be car- 
ried out with the help of the vector identity V ■ (Ax B) = 
B ■ V x A — A ■ V x B, the divergence theorem, and 
VxB = p n KS(y). The surface integral at infinity van- 
ishes, and L' is then obtained from 



A z {x,0)K z (x)dx/I 2 



(31) 



where A z (x, 0) and K z (x) are given by Eqs. (25) and 
(27). The current I z is obtained from the integral of 



6 



K z (x) over the center conductor, 

h = f K z (x)dx = —f^a/A), (32) 

J-a ^MO 



where 



h(u)= / f(u')du'= / 
Jo Jo 



(1 



t(f + 1) 



(33) 



The weak dependence of L' m upon A/a, shown by the 
labeled solid curve in Fig. 6, can be understood via Eq. 
(8). Although we are dealing here with a CPW with di- 
mensions b/a = 1, the penetration of magnetic flux to a 
distance of order A to the sides of the narrow slits leads 
to an increase in the effective width of the slit or, cquiv- 
alently, to an increase in the effective value of the ratio 
b/a. Denoting this effective value as (b/a) e ff — l/fc e // = 
(b/a)(l + cA/a), reveals that c is a slowly varying func- 
tion of A/a of order unity. For example, for the case 
of b/a = 1, substituting (6/a) e // = (1 + cA/a) in place 
of b/a in Eq. (8) and equating the result to L' m calcu- 
lated from Eq. (31) yields c vs A/a as follows: (c, A/a) 
= (0.72, 0.001), (0.68, 0.01), (0.58, 0.1), (0.37, 1), (0.20, 
10), (0.14, 100), and (0.12, 1000). 



2. Kinetic inductance 

The kinetic inductance per unit length L' k is defined 
by Eq. (11) and for thin films is 



L ' k= w/l^ 2(X)d " 



(34) 



Its behavior vs A/a is shown by the labeled solid curve 
for L' k in Fig. 6. 



3. Total inductance 



As can be seen from Eqs. (28), (29), (28), (31), (32), 



and (34), the total inductance per unit length, L' 
L' k is given by 



L' = 



7T^ 



4/i(2a/A)' 



(35) 



The solid curves in Fig. 6 exhibit how L' m , L' k , and 
L' depend upon A/a. When A/a 3> 1, the inductance 
per unit length becomes dominated by the kinetic induc- 
tance contribution p A/4a arising from the nearly uni- 
form sheet-current density I z /2a in the center conductor. 



IV. THIN-FILM CPW WITH WIDER SLOTS 
A. General properties 

Our goal is to calculate the inductances of a CPW for 
which the center conductor carries a current I z e lojt , but 
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FIG. 6: The solid curves show the exact inductances per unit 
length (in units of /io) calculated as in Sec. IIIB: L' [total, 
Eq. (35)], L' k [kinetic, Eq. (34)], and L' m [geometric, Eq. (31)] 
for the limiting case of a thin-film CPW with two narrow 
slits (6 — > a). The dashed curves show the approximations to 
the corresponding inductances calculated as described in Sec. 
IVB5. 



as in Sees. II and III, to avoid including the factor e tuJt for 
all quantities proportional to the dynamic current, it is 
convenient to consider the behavior at a time t for which 

e iujt = L 

Consider now a thin superconducting film in the xz 
plane with two slots parallel to the z axis, one at a < 
x < b containing N' positive flux quanta per unit length, 
and one at — b < x < —a containing an equal number 
of negative flux quanta per unit length. The integral of 
the sheet-current density K z (x) over the width 2a of the 
center conductor is I z , and the integral of K z from b to 
oo (or from — oo to —6) is —I z /2. Of course, K z is zero 
in the slots. Once K z (x) is known, the vector potential 
can be obtained from 



A z (x,y) = - 



Mo 
2tt 



In ^{x-x'Y + y 2 K z {x')dx'. 



(36) 

Because K z (x) is also related to A z (x, 0) via the thin-film 
London equation, we have 

K z (x) = (2//i A)[$'- A z (x,0)}, \x\<a, (37) 
= -(2/fi A)A z (x,0), \x\>b, (38) 



The magnetic flux-density compo- 



where $' = N'< 
nents can be obtained from B x (x,y) = dA z (x,y)/dy 
and B y (x,y) = —dA z (x,y)/dx. These fields have the 
properties that B y (x, 0) is an odd function of x that di- 
verges at x = ±a and x — ±b but decreases to zero 
as | re | — > oo. B x {x,y) is an even function of x but an 
odd function of y with a discontinuity across the film at 
y = 0, given by B x (x, 0+) - B x (x,0~) = -/j, K z (x). For 
p = \J x 1 + y 1 3> A and p 3> 6, screening by the super- 
conducting film causes the magnetic flux density at large 
distances to look like that produced by a dipole of mag- 
netic line charges. In the limit as A/a^ 0, all the fields 
reduce to those discussed in Sec. II. 
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The geometric inductance per unit length L' m is defined 
via Eq. (30) and can be obtained either via Eq. (31) or 
the double integral 



T' — 



Mo 

271-/2 



/OO pC 
-oo J — ( 



In \x - x'\K z (x)K z {x')dxdx' . 



(39) 

The kinetic inductance per unit length L' k is defined by 
Eq. (11) and can be obtained as in Eq. (34). As can 
be seen from Eqs. (31), (34), (37), and (38), the total 
inductance per unit length, L' = L' m + L' k is then given 
by V = &/I z . 



B. Approximate solutions for a CPW 

Numerical solutions of the thin-film London equation 
and Maxwell's equations would be required to calculate 
accurately the fields and inductances for finite gap widths 
b — a with b > a. However, it is possible to obtain rea- 
sonably good approximations to the current and field dis- 
tributions for finite A using the following expressions for 
the sheet-current density, 



K z (x) 



2A 



y/ (a 2 — p 2 x 2 )(b 2 — p 2 x 2 ) 
= 0, a < \x\ < b, 
2A 



, \x\ < a, (40) 
(41) 

, |z| > b, (42) 



y/ (x 2 — p 2 a 2 )(x 2 — p 2 b 2 ) 

where the dependence of the parameter p (0 < p < 1) 
upon A/a is determined as described below. This ap- 
proximation automatically satisfies the requirement that 
J°° K z [x)dx = for all p. When A/a — > 0, we have 
p — >• 1, and K z (x) is given by Eqs. (3)- (4). The constant 
A is given for any p by 



A = 



bpl z 



4F(sin _1 p, a /b\ 



(43) 



where F((j>, k) is the elliptic integral of the first kind of 
modulus k. When p — > 1, F(n/2,k) = K(k), and we 
recover Eq. (5). When k — > 1, F(n/2,1) — tanh _1 p. 
When A -> oo, p -> 0, F(sin _1 p, k) -> p, and A -> bI z /A. 



1.000 
0.500 




1000 



FIG. 7: The parameter p and 1 — p vs A/a (solid curves) 
and the approximations of Eq. (44) (dashed) and Eq. (45) 
(dot-dashed). 




FIG. 8: L' m /no, the geometric inductance per unit length in 
units of fj,o, vs b/a, calculated from Eq. (39) and Eqs. (40)- 
(43) for A/a = (solid, p = 1), 0.1 (long dash, p = 0.932), 
0.3 (dashed, p = 0.817), 1 (dot-dashed, p = 0.587), and 10 
(dotted, p = 0.205). 



It is assumed here for simplicity that the dependence of 
p upon A/a for b > a remains the same as for b = a. 



2. Geometric inductance 



1. p vs A/a 

The dependence of p upon A/a shown in Fig. 7 is here 
obtained by equating the average over the center conduc- 
tor (—a < x < a) of [K z (x) — I z /2a] 2 calculated from Eq. 
(40)] to the average of [K z (x) — I z /2a] 2 calculated from 
Eq. (28) for the case of b = a. For large and small values 
of A/a, the numerically calculated values of p and 1 — p 
are well approximated by 

p w 0.63/v/A>, A/a>l, (44) 
l-p w 0.67 A/a, A/a < 1. (45) 

(46) 



Shown in Fig. 8 are numerically calculated values of 
L' m vs 6/a obtained from Eqs. (39)-(43) for various val- 
ues of A/a. For A/a = (p = 1), the results reduce 
to those calculated in Sec. II and shown in Fig. 3. For 
A/ a < 0.01, the curves of L' m vs b/a are practically indis- 
tinguishable from the solid curve, and for A/a > 10, they 
are practically indistinguishable from the dotted curve. 

Note that when b = a, although L' m = for A/a = 0, 
L' m is finite for all A/a > 0, as shown in Fig. 6. This is 
due to penetration of magnetic flux through the film to 
a distance of order A on both sides of the two slits. As 
discussed in Sec. IIIB 1, L' m can be estimated from Eq. 
(8) upon replacing (b/a) — 1/k by (6/a) e // = 1/^e// = 
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FIG. 9: L' k /[io, the kinetic inductance per unit length in units 
of (j,o, vs b/a, calculated from Eq. (47), where p is calculated 
for each A/a as in Sec. IV B 1 (see Fig. 7). 



(b/a)(l + cA/a), where c is a slowly varying function of 
A/ a of order unity. 



3. Kinetic inductance 

The kinetic inductance is obtained by carrying out the 
integrals of K 2 (x) in Eq. (34). Note that with Eqs. (40)- 
(42), these integrals are convergent, so that, in contrast 
to the behavior discussed in Sec. II B, no additional cut- 
off length needs to be introduced. The integrals yield 
the following convenient analytic result for the kinetic 
inductance per unit length: 



L 'k = ^9k P {k,p), 



(47) 



where 



9k P (k,p) = 



(fc+p 2 )tanh 1 p - (I + kp 2 ) t&nh 1 (kp) 



p(l-fc 2 )[tanh -1 p] 2 



and k = a/b. Note that gk P (k, 1 — e) — > g e (k,e) in the 
limit as e — > 0. 

Figure 9 shows L' k vs b/a for A/a ranging from 0.0003 
to 100. 



5. Comparison of approximate and exact solutions for two 
narrow slits 

We now examine the extent to which the inductances 
calculated exactly in Sec. Ill B for two narrow slits can be 
approximated by the inductances calculated using Eqs. 
(40)- (42). Shown as the dashed curves in Fig. 6 are plots 
of the approximate inductance contributions for b = a. 
The kinetic inductance per unit length L' k was obtained 
from Eq. (47) using Eq. (48), which in the limit k = 
a/b — > 1 reduces to 



_ p(l+p 2 )-(l^p 2 ) 2 tanh l p 
9k ^ P) - 2p(l-p 2 )[tanh-V] 2 ' 



(49) 



and p was determined for each value of A/a using the 
procedure described in Sec. IV B 1. The geometric induc- 
tance per unit length L' m was calculated by numerically 
integrating Eq. (31) using A z (x) and K z {x) obtained 
from Eqs. (36) and (40)- (42), and the total inductance 
per unit length was obtained from the sum V = L' rn + L' k . 

As can be seen from Fig. 6, the approximate solutions 
quite accurately reproduce the exact results for b = a 
for A/a < 1. For large A/a, however, the approximate 
solutions deviate from the exact results chiefly because 
the magnitude of K z given by Eq. (42) is too large for \x\ 
slightly larger than a, and it decreases too rapidly with 
|a;|. The exact result for L' k in the limit as A/a — > oo 
arises entirely from the integral of K z (x) 2 over the cen- 
tral conductor, but the approximate calculation in the 
same limit gives contributions to the integral in Eq. (34) 
from \x\ > a equal to 1/3 of the contribution from the 
central conductor. As a result, the approximation to L' k 
in the limit as A/a — » oo is too large by 33%. The de- 
viation of the approximation to K z (x) from the exact 
K z (x) for |x| > a is also responsible for the fact that the 
approximation to L' m is too small, as shown in Fig. 6. 

Examination of the contributions to L' k from the inte- 
grals of K z (x) over \x\ < a and |x| > b reveal that the 
ratio of the contribution from \x\ > b to the contribution 
from \x\ < a decreases as b/a increases. This suggests 
that the accuracy of the approximate solutions for L' , 
L' k , and L' m when b > a is better than what is shown in 



(48) Fig. 6 for the limit b -> a. 



ACCOUNTING FOR NORMAL-FLUID 
RESISTIVE LOSSES 



4- Total inductance 

The total inductance per unit length is given by the 
sum V = L' m + L' k . Comparison of Figs. 8 and 9 reveals 
that L' k <C L' m when A <C a, L' k <~ L' m when A <~ a, and 
L' k 3> L' m when A»a. Thus, for all b/a, L' is dominated 
by L' m when A«n and by L' k when A » a, similar to 
the behavior of the exact solutions for V vs A/a in the 
limit 6 — > a shown by the solid curves in Fig. 6. 



To deal with the resistive losses at frequencies of in- 
terest to electronics applications, it is useful to make 
use of the well-known two-fluid description, 30,34 in which 
the local current density is the sum of two contributions, 
the superfluid current density and the normal-fluid cur- 
rent density j = j s + j n . The supercurrent density is 
j s = —A s /hq\ 2 , where A s is the gauge-invariant vec- 
tor potential (see Sec. Ill) and A is the temperature- 
dependent weak-field London penetration depth. Be- 
cause the local electric field obeys E = —dA s /dt, we 
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have 



E 



(50) 



when all fields vary as e luJt . The normal-current density 
is 



3n 



E. 



(51) 



where <r„/ = l/p n f is the normal- fluid conductivity, 
i.e., the conductivity of the thermally excited quasiparti- 
cles. a n f is strongly temperature-dependent, being neg- 
ligibly small at temperatures T <C T c and reducing to 
the normal-state conductivity cr„ at T = T c . Note that 
j = aE, where a = g\ — ia 2 = a n f — i/p, cj\ 2 is the 
complex conductivity. 30 

Consider the simple case in which the superconductor 
can be approximated by a semi-infinite half-space x > 
and is subjected to a parallel applied ac field in the y 
direction of amplitude Ho. Making use of Ampere's law 
j = V x H and Faraday's law V x E = —p,odH/dt, we 
find that the induced current j — zj z obeys 



d 2 j z 
dx 2 



\ 2 



(52) 



whose solution is j z = — {Ho/\)e~ x 'e . Here the single 
material-dependent length scale in the problem is A, the 
complex penetration depth 35 ' 36 defined (in the absence 
of vortices) by 



1 1 2i 



(53) 



where <5„/ = (2//Uoajcr„/) 1 / 2 is the normal- fluid skin 
depth. Note that S n f is very large when T <C T c 
that it reduces to the normal-state skin depth 5 n = 
{2/^LO(j n ) 1 ' 2 at T = T c . 

Now consider the case at hand, a CPW consisting of a 
film of thickness d < 2|A|, with central conductor \x\ < a 
carrying current I z e lu)t , slots at a < \x\ < b, and ground 
planes carrying the return current equally on both sides. 
Because the current density j z is very nearly uniform 
across the film thickness, it is convenient to examine the 
behavior of the sheet-current density K z — j z d. The 
analog of j = aE is 



K r = — 



2i 



-E z = Y u E z = (G nf -^—)E z , (54) 



/iqujA/ 



where A w = 2\ 2 /d, the two-fluid version of the 2D screen- 
ing length, is defined via 



1 



I i 



(55) 



A n f — S 2 f/d = 2/jj,ouiG n f is the normal- fluid 2D screen- 
ing length, Y u is the complex sheet admittance, G n f = 
a n fd — l/Rnf is the normal-fluid sheet conductance, 



and R n f — p n f/d is the normal- fluid sheet resistance 
or normal- fluid resistance per square. Inverting Eq. (54) , 
we obtain 



E z 



K z = Z^K Z 



(56) 



where Z w is the complex sheet impedance. 

In the limit as T — > 0, when R n f -4- and A n f -4- oo, 
we obtain Y u -4- —2i/poujA, Z u — > i^o^A/2, and the 2D 
screening length becomes the purely real quantity A w — > 
A = 2X 2 /d, discussed in Sees. II-IV. 

In the limit as T — > T c , when A — > oo and the normal- 
fluid quantities reduce to their normal-state counter- 
parts, we see that Y u — > G n , — > R n = p n /d, and 
the 2D screening length becomes the purely imaginary 
quantity A w — > — iA n = — 2iR n / fiQiu. 

Starting from Eq. (54) and making use of Faraday's 
law and the Biot-Savart law for the general case, we find 
that K z in the CPW obeys not a second-order differen- 
tial equation such as Eq. (52), but the following integro- 
diffcrential equation (omitting the time-dependent factor 
e iut m £ e j^ q uan tities), 



A* dK z (x) 
2 dx 



- H y (x,0) 



P 
27 



J — ( 



K z (x') 



x — x 



T dx', 



(57) 



where P denotes the principal value, and the integral is 
over all x' except the slots (a < \x'\ < b), where K z {x') = 
0. K z is subject to the constraint that /" a K z {x)dx = I z . 
The only material-dependent length scale in Eq. (57) is 

Note, however, that K z {x) and H y (x, 0) vary spatially 
over large distances even when A w = 0. In that case, the 
solutions of Eq. (57) are given exactly by Eqs. (2)-(4), 
which show that K z (x) diverges at the boundaries of the 
slots at \x\ = a and b. For finite A w , although dK z (x)/dx 
and H y (x,0) still have logarithmic divergences at \x\ = a 
and 6, K z (±a) and ^(±6) are both finite. When |A W | <C 
a, the solutions for K z (x) and H y (x, 0) deviate from the 
A w = solutions within roughly |A W | of the boundaries 
of the slots at |x| = a and b. 

The resistive losses in the metallic components of a 
CPW are characterized by the attenuation constant 37 
a c = (—dP/dz)/2P, where —dP/dz is the time-averaged 
power dissipated in ohmic losses per unit length and P 
is the time-averaged power flow along the wave guide. 
From Eq. (56) we obtain 

/oo 
\K z (x)\ 2 dx, (58) 
-oo 

and P — I 2 Z /2, where Z is the characteristic 
impedance of the guide. 

The integral in Eq. (58) is essentially the same as that 
needed to calculate L' k in Eq. (16). To evaluate it ap- 
proximately for |A W | <C a, we again cut off the integrals 
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at x = ±(a — 5 a ) and x = ±(6 + 6b), where 5 a and 6b are 
chosen such that the magnitudes of the current in the 
central conductor and the return current in the ground 
plane are equal. As in Sec. II B, this leads to the condi- 
tion that 6 a /a = Sb/b = e « 1, where e ~ |A w |/a. The 
result for the attenuation constant is 



a c = 



Rnf 



8aZ 



1 



In 



(l-k)K 2 (k) Le(l + fc) 



2(1 - k) 



(59) 



where k — a/b. The factor i?„/|A w /A„/| 2 reduces to zero 
in the limit T -» and to R n in the limit T -» T c . 

References 6 and 38 derived similar but more compli- 
cated expressions for a c containing two logarithmic terms 
arising from the unequal-current assumption of a single 
cutoff length (called the stopping distance A) at both 
\x\ = a and b. Despite numerical differences when b > a, 
the equal-current result given in Eq. (59) is in good agree- 
ment with the results in Refs. 6 and 38 in the thin-film 
limit d < 2|A| when (b — a) <C a. 

As noted in Ref. 2, the total microstrip losses ar are 
the sum of both the conductor loss a c and the dielectric 
loss ad (not considered here), and both these contribu- 
tions must be included to determine the quality factors 
of cavities built from CPWs. 



VI. DISCUSSION 

In this paper I analyzed the geometric, kinetic, and to- 
tal inductances and the attenuation constant for a thin- 
film superconducting coplanar waveguide (CPW) res- 
onator. I pointed out several limitations of expressions 
that were derived assuming negligible penetration of per- 
pendicular magnetic self-fields into the film. Using a the- 
oretical approach that guarantees equal currents in the 
central conductor and the ground planes, I derived new, 
simpler expressions for the kinetic inductance and the at- 
tenuation constant valid for small values of the ratio A/a 
or |A w |/a. I derived exact results for the idealized case of 
a CPW in the limit as b — > a for arbitrary values of A/a. 
I then showed how to estimate the inductances using an 
analytic approximation that reduces to the exact results 



for the current density in the central conductor in the 
limits A/a —¥ and A/a — > oo. 

For the analysis of the kinetic inductance in this paper, 
I assumed for simplicity that the CPW carries sufficiently 
low currents that suppression of the superconducting or- 
der parameter is negligible and A is simply the weak-field 
London penetration depth. However, it is well known 
that application of a sufficiently large bias current can 
suppress the order parameter, making the effective pen- 
etration depth, the kinetic inductance, and the electri- 
cal resistance current-dependent, the changes in all these 
quantities initially increasing quadratically with the bias 
current. 39 ' 44 Similar effects can be produced via induced 
currents flowing in response to an applied perpendicu- 
lar field. 45 These properties make it possible to fabricate 
a variety of devices taking advantage of the associated 
nonlinearities. 24 ' 30 ' 39 - 43 ' 45 ' 46 

To treat current-dependent and associated nonlinear 
effects was beyond the scope of this paper. Within 
Ginzburg-Landau theory, the degree of order-parameter 
suppression scales approximately as (j z /jd) 2 , where j z is 
the local current density and jd is the Ginzburg-Landau 
depairing current density. 25,39 ' 44 As seen in Figs. 1 and 5 
and in Eqs. (3)-(4), (27)-(29), and (40)-(42), the current 
density in the CPW generally is strongly dependent upon 
x, which would make calculations of the current depen- 
dence of both the kinetic inductance and the dissipation 
very complicated, except when A/a ^> 1, when the cur- 
rent density in the central conductor becomes practically 
uniform. 
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